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Abstract: ATP delivery and its usage are achieved by cycling of respective intermediates 
through interconnected coupled reactions. At steady state, cycling between coupled 
reactions always occurs at zero resistance of the whole cycle without dissipation of free 
energy. The cross-bridge cycle can also be described by a system of coupled reactions: one 
energising reaction, which energises myosin heads by coupled ATP splitting, and one de- 
energising reaction, which transduces free energy from myosin heads to coupled actin 
movement. The whole cycle of myosin heads via cross-bridge formation and dissociation 
proceeds at zero resistance. Dissipation of free energy from coupled reactions occurs 
whenever the input potential overcomes the counteracting output potential. In addition, 
dissipation is produced by uncoupling. This is brought about by a load dependent 
shortening of the cross-bridge stroke to zero, which allows isometric force generation 
without mechanical power output. The occurrence of maximal efficiency is caused by 
uncoupling. Under coupled conditions, Hill's equation (velocity as a function of load) is 
fulfilled. In addition, force and shortening velocity both depend on [Ca + ]. Muscular 
fatigue is triggered when ATP consumption overcomes ATP delivery. As a result, the 
substrate of the cycle, [MgATP 2- ], is reduced. This leads to a switch off of cycling and 
ATP consumption, so that a recovery of [ATP] is possible. In this way a potentially 
harmful, persistent low energy state of the cell can be avoided. 



Keywords: energetic coupling; zero resistance; cross-bridge cycle; muscular efficiency; 
muscular fatigue 
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1. Introduction 

In cellular metabolism, energy transductions are brought about by coupled reactions. The network 
of energy metabolism is organised in such a way that cycling of respective intermediates, like protons 
in oxidative phosphorylation (OP), or [Pi], [ADP], and [ATP] in the ATP cycle, is ensured. As was 
shown previously [1], entropy production during steady state cycling must be zero. This follows from 
the fact that the line integral taken around a closed path is zero if the integrand is an exact differential. 
This latter constraint is always fulfilled for potential functions like electro-chemical potentials or 
affinities. It will be shown that this steady state cycling between coupled reactions is associated with 
the occurrence of negative resistances. 

In a muscle fiber, mechanical power output is coupled to ATP splitting. How this is achieved is not 
fully understood, although there has been great success in many field endeavours in muscular research 
such as the structure of the contractile apparatus and its functional correlates [2-6]. Since Huxley's 
widely accepted sliding filament theory [7,8], the cross-bridge cycle is of central importance, 
especially in the functional aspects of contraction. This cycle must contain the reactions of free energy 
transduction from chemical (ATP splitting reaction) to mechanical energy (actin movement against a 
load force). From the overall reaction, contractile efficiency can be obtained by relating mechanical 
power output to the dissipation function of ATP splitting, where the mechanical power is given by the 
product of the force exerted by the load and the shortening velocity. Experimental results show [9-12], 
that when efficiency is expressed as a function of v, a curved line with a maximum is obtained. From 
non-equilibrium thermodynamics (NET, [13]), it is well known that uncoupling is necessary to 
generate a maximum in efficiency plots (efficiency against reduced force ratio). Thus, to yield such a 
maximal efficiency, any description of the cross-bridge cycle on a thermodynamic basis must contain an 
uncoupling mechanism, which uncouples the transduction of free energy from ATP splitting to 
actin movement. 

To describe the cross-bridge cycle in terms of the new flux equations published recently [1], the 
cross-bridge cycle has to be formulated in relation to this formalism, which combines the basics of 
NET [13-16] with Michaelis-Menten-like kinetics of enzyme-catalysed reactions [17]. It will be 
shown that Hill's equation describing muscular performance [18,19] can be easily deduced by 
applying the new flux equation. 

When compared with other approaches to the energetics of the cross-bridge cycle, the main 
particularity of the present work may be the fact that this cycle is connected here to energy metabolism 
of the muscle fiber, i.e., to ATP producing and consuming reactions. The generation of mechanical 
energy from the free energy of ATP splitting is treated here as one of the parallel reactions of the 
sarcosol consuming ATP delivered in fast fibers, mainly from glycogenolysis or glycolysis, 
respectively. This integration into the cell's energy metabolism makes it possible to inspect some 
variables like ATP and its reaction products and species at high mechanical power output. In addition, 
concentration changes in metabolites and ions like creatine phosphate, lactate, H , and Mg , are of 
interest under these conditions. This is achieved by formulating, in particular, the ATP splitting 
reaction according to Alberty [20] as a function of both [H ] and [Mg ]. 

It is the aim of this study to elucidate cycling between coupled reactions, and if such cycling is also 
involved with the cross-bridge cycle and force generation. In addition, the consequences of uncoupling 
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on power output and efficiency will be shown. For this, a formulation of the cycle in terms of the 
above mentioned new flux equation had to be derived. 

The phenomenon of muscular fatigue at the cellular level occurs when ATP consumption exceeds 
ATP delivery [21-25]. Under such conditions drastic changes in many metabolite and ion 
concentrations can be expected. The results of simulations will show to what extent these changes may 
contribute to fatigue, and if this phenomenon can be explained by such changes alone. 

2. Results and Discussion 

2.1. How Negative Conductances Are Generated 

In a previous article [1] it was shown that at steady state all affinities and dissipation functions of 
closed pathways associated with coupled in series reactions must vanish. In the following, it will be 
demonstrated that the overall resistance (=1 /conductance) of such cycles must also be zero. As a 
consequence, the existence of negative conductances (or resistances) has to be called for. 

According to [1] a coupled two-flux-system can be described as: 

J r = L c ((A, +l)A 1+ A 2 ), and J 2 = L c (A, + [X 2 + 1) A 2 ) (1) 

J j and J 2 designate fluxes through affinities A } and A 2 , respectively, and L c represents the coupling 
conductance. Under totally coupled conditions (\= X 2 = 0) both fluxes are equal. The dissipation 



function, & , of a coupled process is composed of two parts, 0 1 for the output, and @ 2 for the input 
reaction, with: 

0 = 0 1+ 0 2 (2a ) 
&j =JjAj, and & 2 =J 2 A 2 (2b) 
0j = L c (Aj + A 2 )Aj , and 0 2 = L c (A } + A 2 )A 2 (total coupling). (2c) 

because A i usually is negative, 0 l must also be negative. Expanding the right hand terms yields: 



f A±A^ 



A) , and 0 2 = L c 



r A 1 +A 2 ^ 

V A 2 J 



(2d) 



The term: 



L c i=k 



r A ] +A 2 ^ 



(2e) 



represents that partial conductance of L c , which is associated with A /, while 



L c2 =k 



r A±A 2 } 

A 2 



(2f) 



belongs to A 2 . They relate to the usual different forms of energy being processed through the coupling 
reaction. Obviously, when A } is negative, L cl must also be negative to yield a negative 0 1 . 
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The same result can also be derived by starting from flux equations, yielding: 

A, + A 2 _ 

J 



f A . A \ ( A+A^ 

A 1 , and J 2 = L 



- A 2 



A, (3) 



x 2 



So, to yield a positive Ji, L cI has to be negative for a negative Aj. 

A j and A2 are in series, hence, L c can be regarded as the equivalent conductance of both in series 
conductances L cl and L c2 , yielding: 



L+J_ c 4 1 4 (4) 

L cl L c2 Aj + A 2 Aj + A 2 
These theoretical results are confirmed by simulations. 
2.2. Conductances in Cycles between Coupled Reactions 

In a reaction sequence in which two coupled reactions in series are involved, the output force of the 
first reaction, A\ , may be used by the second reaction as an input force Af ( A 2 is the input affinity, 

Af denotes the load affinity). In such a cycle between two coupled reactions, both forces must be 

equal but of opposite sign. The output power of the first reaction delivers the input power for the 
second reaction by flowing through A[ , and Af = -A[ , and back to A[ (at zero power). At steady 

state, fluxes through A[ and Af are equal, and hence, both dissipation functions of the cycle, <S>\ and 

0" , must vanish. From &{ = -0f , and L' cl (4) 2 = -L n c2 (A"f , L' cl = -L" 2 is obtained. That is, the 

partial conductances of two coupled reactions in series are opposite and equal, if under conditions of 
steady state cycling the magnitudes of the output force of reaction / and the input force of reaction // 
are equal. 

0\ = -0 2 can also be expressed in terms of the steady state flux of cycling and of resistances 

R' cl = \/L' cl , and R" = \/L" 2 , yielding R' cI J 2 + R" 2 J 2 = 0 . It follows that at steady state cycling between 

two coupled reactions, the sum of the resistances in that cycle must vanish, i.e., that the steady state 
flux through a cycle between two coupled reactions always occurs at zero overall resistance. Cycling is 

driven solely by A 2 + Af ( Af negative). Because both reactions are coupled, the conductance 

(resistance) of the whole cycling process is brought about by both partial (in series) conductances 
associated with the input and load affinity, respectively. 

In oxidative phosphorylation (OP), as described in detail in [1], a proton cycle is generated over the 
inner mitochondrial membrane. At steady state, coupled outward proton pumping by redox (NAD re d 
and FAD re d) reactions of the respiratory chain (J NA and Jfa) equals the back flow of a given fraction 
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(Qh) of protons through ATP synthase (Jsy) and ATP/ADP exchange (Jae) plus H/Pi symport (Jp,). 
Both partial conductances, 



Lri -L r 



v R Ajx H 



r v R Afi H +A° R ) {-A ATPc +(-v s+1 A£i H ) 
, diiu i-, SY2 -J^ SY 



(5) 



are opposite and equal. For the first reaction (index R) A{ = v R A/u H , and for the second 

A" = -v s+1 Aju H , with v R = v s+1 (total coupling; A° R = affinity of Jna plus J fa', v r = v s+1 = 4 

protons/extent of reaction). 

The remaining fraction of protons (Qr H ) flows back (driven by -Aju H ) through several parallel 

conductances given by the proton leak flux, J PL , mitochondrial Na + /Ca 2+ exchange (with Na + /H + 
contracted to H + /Ca 2+ exchange), Uhce, and the malate-aspartate shuttle, Jms- The partial conductance 
of this residual proton efflux and the sum of conductances of back flowing fluxes, are also of opposite 
equality. 

Analogously, partial conductances of ATP cycling through the potentials of mitochondrial ATP 
(ATP m ) production plus ATP transport (contracted to A\ ), and of cytosolic ATP splitting ( A" ) can be 

formulated. Opposite equality of partial conductances is also fulfilled for this cycle (the above results 
were obtained by using the simulation SIMgiox from reference [1]). 

For a further illustration, an analytically solvable example is given in the Appendix section. Simple 
electric circuits consisting of one battery connected to an outer conductance, or of two batteries in 
series, are analysed. These examples show very clearly the behavior of coupled in series reactions. 

Further evidence of such an equality of conductances comes from the known fact that for a coupled 
reaction with an attached load, conductance matching {L L d = L c ) is needed to achieve a maximal power 
output [1]. At total coupling, the output power is given by: 

®i=P mtt =LM 1 +A 2 )A 1 (6) 

The maximal P out is found by differentiation with respect to the variable A\, while A 2 remains 
constant, and by setting the derivative equal to zero: 

dP 

-^L = L c (2A ] + A 2 ) (6a) 
dP 1 

"^ = 0 leads to A lmax =--A 2 (6b) 

Inserting A ]max into equation (2e) yields L cl = -L c and because -L cl = L Ld , it follows L Ld = L c . 

In addition to cycling at the inner mitochondrial membrane, other types of cycles occur in 
metabolism. Especially in skeletal muscle cells, the phosphofructokinase (PFK) reaction in conjunction 
with the fructose- 1 ,6-biphosphatase (FBPase) operating anti-parallel represent a substrate cycle, which 
may control the pathway of glycolysis (GLY) more sensitively than would be possible by PFK alone. 
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In this cycle, fructose- 1,6-biphosphate (FBP), which is produced by ATP-coupled formation from 
fructose-6-phosphate (F6P), is cycled back via FBPase to F6P. However, usually both fluxes are not 
equal. Also to demonstrate the opposite equality of partial conductances for this kind of cycle, only 
equal fluxes can be used for this purpose. 

As a further example, the phosphocreatine shuttle will be considered. The creatine kinase (CK) 
reaction can also be regarded as a coupled reaction. Here, ATP splitting powers phosphocreatine (PCr) 
formation from creatine (Cr), which may proceed near equilibrium. As described in detail in reference [1], 
ATP is shuttled between locations of ATP formation (for instance in the inter-membrane space in 
mitochondria) and locations of high ATP demand like myofibrils. By analogy to an electric circuit 
built by two in series batteries with an outer circuit conductance (see Appendix (A4)), the output 

affinity of PCr formation in the inter-membrane space of mitochondria corresponds to A\ with 

associated L' cl , whereas the affinity of the reverse reaction in myofibrils corresponds to A" (with L u c2 ). 

To ensure diffusional flow of PCr and Cr between both locations, an additional driving force 
(corresponding to U e ; see (A4)) with associated conductance must be present. Under such conditions 
partial conductances do not match. Only when the additional conductance corresponding to the 

diffusional process (L e ) is added to l! cl does this sum become opposite and equal to L u c2 , as is shown in 

(A4). L e depends greatly on structural features. So, to achieve a high diffusional conductance, 
diffusional paths must be as short as possible, which in turn requires a high grade of structural 
organization [26-28]. 

It seems worth mentioning that coupled systems like pump and leak cycles are often not in a steady 

2_|_ 

state. For instance, steady state cycling through sarco/endoplasmatic reticulum Ca ATPases 
(SERCA) and Ca release channels of the sarcoplasmatic reticulum (SR) breaks off during activation 
of contraction. There is an enormous Ca efflux through release channels; meanwhile the pumping 
rate of SERCAs may be low. Under these conditions, respective conductances may greatly differ; 
however, when a new steady state cycling is reached, the partial conductance of SERCA must be 
opposite and equal to the conductance of the Ca release channels. The opposite has to be expected, 
when release channels close again, and the Ca 2+ pumping rate exceeds the release rate. 

2.3. From Chemical Potentials to Mechanical Force Generation 

In striated muscle cells like ventricular muscle cells (VMs) or skeletal muscle fibers (SMFs), force 
generation as well as shortening is brought about by the cyclic action of cross bridges. It is a known 
fact that this process is powered by ATP splitting. The underlying mechanism of the energy 
transduction process, however, is not completely understood. Here, a thermodynamic description of 
the cycle is derived using a formalism recently published [1]. It takes into account the basic energetics 
of enzyme-catalysed reactions, which states that the overall affinity of the catalysed and non-catalysed 
processes must be equal. For an enzyme-catalysed reaction like: 

S + E^SE^SP^S + P 
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(S = substrate, E = enzyme, ES = enzyme-substrate complex, P = product), this means that at steady 
state the sum of the affinities of substrate binding, transition, and product release must yield the 
affinity of the non-catalysed reaction, which is given by the reaction affinity of all involved 
compounds in the bulk solution: 



A r = Kiln 



0 



\ f 

+ RT/n 



J 



v [ES]c 

or, after contraction of the first two terms 



[ES] 
[Epf 



+ RT/n 



J 



o \ 



[P][E] « 



(7a) 



A r = Kiln 



[EP] 



yielding, 



A = RT/n 



+ RT/n 



[s] c 



[EP]i ... 
[P][E] « 



(7b) 



[P] 



(7c) 



With K' r = K' B x K' T x K' R (K' B , K' T , and are equilibrium constants of the binding, transition, and 

release reaction, respectively, whereas K' r denotes that of the non-catalysed reaction). 

An analogous reaction sequence is used here to describe the cross-bridge cycle. The following cycle 
is given in chemical notation, i.e., the charges of involved species are taken into account. The cycle 
begins with the splitting reaction of the de-energised actomyosin complex (A-M) by MgATP " in the 
diffusional space of myofibrils: 

A- M+ MgATP 2 ^ A 2 +M 2+ » MgATP 2 " Rj 

This first reaction yields dissociated actomyosin with MgATP 2 " bound to myosin (the bold point 
denotes binding to myosin). Two negative charges develop on the dissociated actin, which are 

neutralised by potassium ions, K + , stemming from free MgATP 2- , which is now bound to myosin 

heads. On the dissociated myosin heads, it neutralises both emerging positive charges. This first 
actomyosin dissociation and binding of MgATP 2- to myosin is followed by ATP splitting on the 
myosin heads. This transition reaction is described by 



M «MgATP ^ M .MgADP .H 2 PO 



R2 



It is coupled to the formation of energised myosin ( M 2+ ), which is characterised by a tilting of the 

myosin head from a more bent arms position by an angle of about 60° towards the respective Z disc, so 

that now the myosin head builds a right angle with the opposing actin filament. M 2+ contains free 

energy from reaction R2 as conformational energy. 

The force generating stroke of the myosin head is triggered by the association reaction to form the 
energised actomyosin complex (cross-bridge): 
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A 2 ~+M 2+ .MgADP~.H 2 P0 4 ~ ^{A-M.MgADP~.H 2 P(V} 2 ~ R 3 

Because of uncompensated charges, the resulting intermediate in curly brackets lacks firmness. A 
more stable conformation is obtained by the last reaction of the cycle, which restores the electro- 
neutrality of the cross-bridges: 

{A-M.MgADP~.H 2 P0 4 ~} 2 ^A-M + MgADP +H 2 P0 4 R 4 

The existence of a less stable interaction of myosin with actin has been shown previously [29-31]. 
In a coupled reaction, de-energised actomyosin is restored by dissociating first MgADP and 

H 2 P0 4 from cross-bridges and then by releasing the stored conformational energy. During this 

reaction the cross-bridge tilts back by 60° towards the sarcomere centre, whereby free energy is 
transferred to the actin filament as mechanical energy. 

From the above scheme (Ri to R4) two fluxes can be obtained, which are responsible on the one 

hand for the production of dissociated and energised myosin heads (J En ), and on the other hand for the 

formation of cross-bridges and subsequent mechanical force generation by stroking (J Str ). At steady 

state, a certain fraction of myosin heads of a half-sarcomere exists in a dissociated and energised state 
M 2+ «MgADP «H 2 P0 4 ([MH En ]) , while the residual fraction interacts as cross-bridges with actin. The 

resulting fluxes are given by: 



4 n = RT In 



K B 



ref 



and 



(8b) 



A" Str = RT In 



^En=L E n(4 d n+4n),^ (8a) 

[MgATP 2 -] 
[MH En ] 

J Str = L St r ( A st + Kr ) > With (8c) 

[MHJ[C " ] - .(-4!) (Be) 



J 



[MgADP"][H 2 P0 4 -] R 

(For a more complete description and definition of reference constants (K ref ) see (A5); complete 
conductances (L En and L Str , respectively) are given in (A14) and (A15).) 

If the constraint *K re ^ = K re ^ is fulfilled, contraction of the affinities of both fluxes yields the 

B R ATP 

overall affinity as required ( A^f r + A ATP ). Here A L E d n (stored as conformational energy) denotes the 
affinity coupled to binding of MgATP 2- to myosin heads ( A P En ), and A L s d tr the affinity which is coupled 
to the power stroke potential ( A p Str ). A L s d tr represents the mechanical work per mole of cross-bridges 
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which has to be overcome during stroking. The quantity J Str x A Str is directly related to mechanical 

power output P Str = F u x v {F L d = load force in Newton (N) of all stroking cross-bridges of a given 

cross sectional area; v = velocity of shortening in m/s related to a given fiber length), which as such is 
conveyed to the surroundings. 

In the present model of the cross-bridge cycle, A ATP is used at two mechanistically and temporally 

separated steps. They are given on the one hand by binding of MgATP 2- and on the other hand by the 
release of MgADP" and H 2 PO/f . Here, most of the free energy of ATP splitting is associated with A P En , 

which by the coupling process on myosin heads is transformed first into A L E d n , and then is delivered as 

A P Str to the power stroke after cross-bridge formation. Therefore, the stroke potential in mechanical 

units (cross-bridge force x stroke length x Na, Na = Avogadro's number) must be equal to A P Str (see below). 

Because ionic species are involved, the reaction sequence of the cycle should be markedly enforced 
by electrostatic interactions. So MgATP binding can proceed only if actomyosin dissociates, whereas 
release of products becomes possible only when at the same time cross-bridge formation occurs. 
Moreover, the conformational change in the myosin head forces it into a new position, which favours 
an interaction with actin at a new actin binding site displaced a certain distance towards the Z-disc. 
During stroking, binding of a new MgATP molecule and detaching of cross-bridges may 
preferentially occur at the end of the power stroke, when cross-bridges form an angle of about 60° with 
the actin filament (see below for uncoupling by stroke shortening). 

The contractile performance of whole muscle and of SMFs is exceptionally well reproduced by 

Hill's equation [19]. This equation relates the shortening velocity v to the mechanical load force F Ld 
which has to be overcome during shortening. 



The above function represents a hyperbola, which fits remarkably well with experimental data 
obtained under isotonic conditions. 

To obtain an equivalent expression from the flux equation J Str , the flux given in mM/s has to be 

converted into velocity with units of m/s. This is achieved by calculating the stroke frequency for a 
given concentration of stroking cross-bridges ([CB] = [CB] to t - [MH En ], in mM) and by multiplying 

with the stroke length l Str (in m) and the number of in series half-sarcomeres Nh s . The result is: 



v 




a + F, 



(9a) 



Ld 




(9b) 
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The above expression describes the shortening velocity as a function of A$f r at constant A P Str . It 

represents a straight line (Figure 1A). Introducing a Michaelis-Menten like inhibition factor associated 
with L Str yields the desired hyperbolic dependency: 



[CB]Km u +A£ 



£m„ LJVj x ^" + ^; ; 



[CB] 



(9c) 
(9d) 



Comparing equations 9a and 9d shows that the constant b of Hill's equation is given by: 



b 



-'Str 



[CB] 



(9e) 



As required, the quotient by which b is multiplied is dimensionless. To yield the shortening velocity 
as a function of force, v(F u ), affinities and Km Ld (both in J/mol) have to be converted into units of 

force. This is achieved by dividing by l Str and by multiplying by the molar number of cross-bridges. 



Figure 1. Flux as a function of load potential at 10.8 uM [Ca ]. A: (grey dots) according 
to equation 9b; (light grey dots) according to equation 9c or 9d; (red line) according to 
equation 11a; (green line) according to equation lib. B: (light grey dots) according to 
equation 9c or 9d; (red squares) results from simulation SlMoLYgen- 




6x10 




6x10 



-Aload [J/mol] 



-Aload [J/mol] 



A\* being negative, F Ld must also be < 0. Expressing shortening velocity as a function of a 
positive variable yields with F Ld = -F^ d 



v(F u ) = -bx 



F P ~F U 
-Kfm Ld+ F + Ld 



(10a) 



Setting -Kfm Ld = a , and -b = b + , gives (Figure 2.) 
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v{F u+ ) = b^ F "- f ''- 



a + F, 



(10b) 



Ld+ 



2+ 



Figure 2. Shortening velocity as a function of load force at two different Ca 
concentrations A:[Ca 2+ ] = 1.08 uM; (light grey dots) according to equation 10b; (red line) 
equation 10b plus uncoupling; (red circles) results from SlMaiYgen versus load force; (blue 
squares) results from SlMciYgen versus load force as sensed by cross-bridges; B: as in A, 
but at [Ca 2+ ] = 0.36uM. 



2x10 



1x10 " 




2- 



8x10 



2x10" 



1x10 




1.2x10 



The latter equation formally represents Hill's equation. In that equation Fq denotes the maximal 
force obtained under isometric conditions, whereas Fp in the latter equation is obtained from the input 

affinity ( A P Str ) of J Str by converting it into units of force (see below for a derivation of F p = F 0 ). 

When efficiency is represented as a function of shortening velocity, the experimental data follow a 
curved line with a maximal efficiency at 0.18 x v max [12]. From NET it is known that such a maximum 

is produced by uncoupling. To create such a maximal efficiency, uncoupling terms have to be 
incorporated into J Str . Variable, load dependent A values ( A(Af£) instead of constant A's ) are 

defined, to preserve the hyperbolic nature of the function. In this way, uncoupling becomes operative 
only when A^f r exceeds a certain value (Figure 1 A). Both flux equations are given by: 



■#=M(^(4S) + l)^+4r) 

Jl-L S ,{A^ + {XI{A^) + \)A P S J) 



(11a) 
(lib) 



These latter equations (for a complete description, especially of conductances, see (A 15)) appear in 
simulations. A values are given as functions of A^f r , e.g., 



^Strf 



1 + exp 



( Aid 

A Str 



A 



f \ 



Id 



'Ld 



1 + exp 



f iLd ~ls 

A Str A Ld 



'Ld 



(12) 
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(see below for a mechanistic interpretation of uncoupling and A values). Conversion to mechanical 
units can then be done in the same way as shown above. 

In Figure 1 the effects of uncoupling are shown. At a load of about — 3.0><10 4 J/mol, deviations from the 
hyperbolic (coupled) curve begin to arise. From the plots it can be seen that uncoupling leads to a shift of 

the intersection with the abscissa to less negative values of A^f r , whereas J p Str is still maintained, even at 

Aj£ = -A p Str , where the coupled flux must be zero and only uncoupled fluxes are possible. 

In the following, an attempt has been made to interpret the above results, which were gained from a 
phenomenological approach, mechanistically by relating coupled and uncoupled fluxes to possible 
cross-bridge actions. 

At Agf r = -A P Str , coupled reactions with associated actin filament movement come to a halt, 

because the driving force has vanished. As already mentioned above, now only uncoupled fluxes can 
occur. Such a situation may also be realised with isometric contraction, which is known to be 
associated with ATP splitting and heat production, but without power output. That is, a mechanism has 
to be found which explains the identity of the isometric force Fq with F P , which was merely formally 

derived from the input affinity A P r by a conversion factor. This is achieved by defining the uncoupling 
mechanism by a shortening of the stroke length l Str of the power stroke. Total uncoupling is reached 
when l Str = 0 . This may be realised under isometric conditions. Free energy corresponding to 

A p r ~ A ATP is delivered to actin filaments as mechanical work, i. e., F 0 xl Str xN A =A p r . Shortening 

may be brought about through splitting of actomyosin bonds before the whole stroke length is 
transferred to an actin filament. When A L s d tr = -A P r , as is realised under isometric conditions, 

actomyosin splitting already occurs at zero stroke length, so that no energy can be delivered to the 
actin filaments. Only force development by cross-bridges during the time interval between bond 
formation and bond splitting is possible under these conditions. This may be achieved by the torque 
every myosin head exerts on an actin filament after bond formation and release of F^PO/f and 
MgADP". The associated force then acts on these filaments, but without being able to bring about 
filament movement, since this is hindered by the equal and opposite load force. Therefore, if the force 
remains constant over the whole stroke length, then F P is equal to the isometric force F 0 . 

When compared with a cross-bridge cycle during contractions, the cross-bridge cycle under 
isometric conditions becomes altered insofar as coupled stroking is impossible; the power stroke 

occurs completely uncoupled, so that all free energy associated with A p r becomes dissipated as heat. 

Moreover, dissipative stroking under these conditions may occur in the presence of bound MgATP . 
The following derivation shows how stroke shortening may be involved with uncoupling. 
Stroke shortening is given by: 
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Al = l 'str -/^ (negative) 



(13a) 



this leads to, 



zl^ r = F Q Al = A P r - A F Str (at constant force) 
Under totally coupled conditions, the input flux is given by: 

J Str ~ ^Str \ A Str + A Str ) 



(13b) 



(13c) 



Uncoupling by stroke shortening dissipates free energy, which can be expressed by a leak 
dissipation function: 



0L=L P StrL (AAif=L 



p 

StrL 



V A Str J 



(13d) 



The leak conductance L P StrL can be replaced by L Str , because this latter conductance may depend 

mainly on the formation mechanism of the actomyosin bond. The stroke reaction associated with 
conformational changes of the myosin head is assumed to proceed at a high conductance, since the 

energising reaction ( J En ), which is coupled to the same conformational change in the reverse direction, 

also proceeds at a very high conductance. So an increase by stroke shortening of a high conductance 
(stroking) in series with a low conductance (bond formation) may be negligible, so that 0 P trL can be 



expressed as: 



( aaP \ 



® StrL =L Str 



AA 



Str 



A p 



K) 2 



Comparing this latter equation with that used in the simulation, 



< L =Lstr4rKr){<)\ yMds: 



(13e) 



( A AP \ 



AA 



Str 



V ^Str J 

The input flux then is given by: 



V htr J 



= 4r(<) 



J Str ~ L Slr iy^str + Ktr ) + Asr 



V A Str J 



'Str 



(13f) 



(13g) 



The output flux is reduced by stroke shortening as if it were uncoupled. The same dissipation 



function L Str (AA P r ) is associated with output reactions, yielding: 



0^ L =L Str (AA P r f=L s 



f A4tr 



V Str 



(13h) 
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and the output flux: 



jLd _ r ( A Ld jP \ j 
J Str ~ ^Str \ A Str + A Str ) + ^Str 



It follows, 



t^ 2 



A Ld 



( a p V 

4 Ld 
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V A Str J 



Aid 

1 Sli- 



V htr J 



(130 



(13j) 



For Al = 0, identical coupled fluxes arise, and for A L * r = -A P Str , both X values are equal. Moreover, 
if A,™ = A P r = 1 , (Al/l Str ) is also equal to 1.0, which means that now isometric conditions do exist. 
From equation (13i) it can be taken that uncoupling by stroke shortening reduces J L s d tr as if there were 
a leak flux through A L s d tr . On the other hand, J P r increases (equation (13 g)) as if there were an 

additional leak flux through A p Str . The above derivations demonstrate that stroke shortening obviously 

leads to the same effects as uncoupling by leak fluxes. It seems justified, therefore, to also describe 
uncoupling by stroke shortening by lambda values, as was done previously mainly in the context of 
oxidative phosphorylation. 

The degree of coupling is given by, 

1 

' ' (14a) 



7(4+1)0^+1) 



[i] 



with above results this yields: 

<lstr 



f < ' A p V 

V ^Str J 



(14b) 



Under the limiting conditions of isometric contraction ( A L s d tr =-A p r ; (zl/// Sr ) 2 = 1.0), q Str is given by: 

1 



^«) + l 



= 0.5 



(14c) 



At loads «-3xl0 4 < A L s d tr < ~ A p r , X p tr (A^f r ) will be smaller than 1.0 but not zero, so that the 

process in this region would proceed at a higher degree of coupling (e. g. at A$ d = -4.5 10 4 , q Str = 

0.833). At values of A L s d r > «-3xl0 4 the process is totally coupled ((Al/l Str ) 2 = 0), that is, cross- 
bridges work at full stroke length. Only this part of the performance curve (Figures 1 and 2) is 
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hyperbolic and fulfils Hill's formalism. Between the intersection {A L s d tr = -4.756x10 ) and A Ld r = -A Str 



, Jgf r formally could be negative, which would mean that actin filaments were moving in the direction 

of stretching. This is, however, impossible, because actomyosin bonds would have to be broken by a 
load force, which is smaller than F 0 . Therefore, in this region of loads, J L s d tr cannot be negative; it must 

remain zero. 

2.4. Power Output and Efficiency 

In experiments, mechanical power output is often represented in relation to shortening velocity. In 

2_|_ 

Figure 3, power and efficiency plots at two different [Ca ]s (1.08 and 0.34 uM, respectively) are 
shown. Respective curves have similar shapes; however, F 0 and v max , and therefore power output 
values, are markedly increased at high [Ca ]. 

Efficiency curves at both [Ca + ]s are nearly identical (Figure 3D). In panel B, efficiency of a totally 
coupled cross-bridge cycle is shown. Under these conditions the curve has no maximum. 

Partial conductances can be calculated from L En , A L E d n , and A P En , as well as from L Slr , A L s d r , and A P Sir 
. All results derived in the above sections could be verified by the simulation (SIM Gi y ge „). So, 

L r „ = — — - — : — , and L^„ = 



'* J_ + J_' Str J_ + J_ (15) 
^Eni L En2 L Strl L Str2 



also, L Enl = -L Str2 is fulfilled, and therefore, cross-bridge cycling at zero resistance. 
In addition, the equality of 

1 1 



^ + J- J_ + J_' (16) 



L En 2 L Str] L En L Str 

which describes the conductance of the whole cycle including coupled inputs and outputs is nearly 
exactly obeyed. 

The overall efficiency of the cross-bridge cycle is obeyed: 

jLd jLd jLd jLd jLd 

n _ A Str _ En A En J Str A Str 

lov j j P .P j P A P (17) 

/4 TP u En ^En J Str ^Str 

as is the overall dissipation function given by: 

0 O v = J En K + A ATP) = < + K + < + K (18) 

2_|_ 

Figure 3D shows efficiency curves at 1.08 and 0.36 uM [Ca ]. They are very similar; their 
maximum lies at about 0.18 v max . Because the appearance of the maximum is caused by uncoupling, 

the coordinates of r] max are highly dependent on uncoupling parameters. 
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Figure 3. Power output and efficiency at two different Ca concentrations. (A) and (C) 
[Ca 2+ ] = 1.06 uM; (B) [Ca 2+ ] = 0.36 uM; C: under totally coupled conditions; (D) (red 

2+ 2~i~ 

squares) efficiency at 1.06 uM [Ca ],(blue circles) efficiency at 0.36 uM [Ca ]. All plots 
are results from SlMoLYge- 
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Calcium Ions and Force Development 

In the previous section it was shown, how shortening velocity depends on Agf r at a given [Ca 2+ ]. 

On the one hand, the driving force is changed by the load potential (see Figure 1, linear dependence), 
and on the other hand the conductance L str depends on A^' r through the hyperbolic inhibition factor. 

At a given [Ca ] s both effects are responsible for the characteristic appearance of the performance 
curve under coupled conditions. 

In the present model of the cross-bridge cycle, interference of [Ca 2+ ] with A p En as well as with L Str is 

necessary. In the latter case, [Ca ] can activate J str through a sigmoid activation factor (A15). This 

takes into account the fact that Ca 2+ binding to troponin C removes the inhibition of cross-bridge 
cycling, so that binding of myosin heads to actin binding sites becomes possible [32,33]. On the other 
hand, [Ca ] is known to strongly activate force development. Here it is assumed that this may be 
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caused by an increase in cross-bridge concentration [CB] . By introducing a [Ca 2+ ] dependent K 
(see (A14)), a sigmoid variation in both [CB] and force F by [Ca 2+ ] can be obtained (Figure 4.). 
Figure 4. Developed force and cross-bridge concentration [CB] and their dependence 

2+ 

on [Ca ]. (red squares) force; (blue circles) [CB]. Notice that at the given dimensioning of 
the right ordinate a matching of results is obtained. 
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At steady state, a certain [CB] is produced by [Ca 2+ ] activated J str , and in addition by [Ca 2+ ] 
inhibited J En (see (A 14)). The inhibition of J En by [Ca + ] is brought about by a decrease of A En with 

2_j_ 

increasing [Ca ]. This is possible because this reaction proceeds at a very high conductance and 
therefore, is close to equilibrium. So already a small variation of the driving force can produce a large 

change in the reaction velocity. In this way, a sensible, [Ca 2+ ] dependent adjustment of [CB] and force can 

be achieved. An elevation of [Ca 2+ ] thus increases both shortening velocity as well as force development. 
The total myosin head concentration ([MHeJ + [CB]) of a half-sarcomere amounts to 656 uM (see 

Methods). At a saturating [Ca ] of 1.08 uM, fluxes J En and J Str are so adjusted as to yield a 
concentration of [CB] = 0.25 ([MH En ] + [CB]), i.e., at this [Ca + ], 25% of myosin heads form cross- 

2"+" 

bridges and thus are involved with cycling and force generation. At [Ca ] = 0.36 uM, only about 3% 
of cross-bridges are engaged, and at 0.09uM [Ca + ], [CB] is further markedly reduced, which means 
that now near resting conditions are reached. 

It seems plausible to suggest that during shortening it is not always the same group of cross-bridges 
that is active, but that, e.g., at 1.08 uM [Ca ], four different groups may alternately be involved with 
contraction. The cycling frequency of an individual cross-bridge would then be much lower than the 
frequency of ATP splitting, which might be advantageous, especially at high velocities. Furthermore, 
an alternating involvement of groups may be absolutely necessary for a smooth shortening. How this 
might be accomplished is so far not known. An involvement of special filaments of the sarcomere 
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cytoskeleton [34-36], which may be responsible for a subtle sensing of load forces and an undisturbed 
takeover of a given load by a new fraction of cross-bridges during synchronous stroking, 
seems indispensable. 

The values of maximal tension (=force/unit area in N/m = Pa, Pa = Pasqual) obtained from 
SlMoiYgen (A16) in the present study are comparable to experimental values. For instance, a value of 
372 kPa (from F 0 = 7.756 x 10" 4 N, [Ca 2+ ] = 1.06 uM, 37°C) found here, seems to be in reasonable 
agreement with about 320 kPa resulting from measurements with a fast-twitch mouse fiber at 25°C [37]. 

The extrapolated value of maximal shortening velocity, vff T = 1.95 um x HS -1 x s -1 ([Ca 2+ ]= 1.06 

uM, A L = 0 J/mol, HS = half-sarcomere) compares to 1.6 um x HS" 1 x s _1 of frog fibers at 0°C [12]. A 
value of rj max of about 50% at about 0.18 v max ([Ca ] = 1.06uM) results from adjustment. It compares 

to the experimental values of 35^45 % for the same value of v for frog muscles at 0 °C [12]. 

All these parameters of contractile performance may, however, be reduced to a certain extent by 
dissipative frictional processes associated with v which are not addressed in the present simulation. 
Such dissipation during fiber shortening may be produced mainly by viscous deformations of 
membranes and the filament lattice. 

2.6. [H*], [Mg 2+ ], and Fatigue 

Enzyme-catalysed ATP splitting by myosin heads is formulated here with respect to the ATP 
species MgATP 2- . By using a reference constant and binding polynomials, an [H + ] and [Mg 2+ ] 

dependent K' ATP of this reaction can be formulated (see A6 and A7). In simulations of fatigue, in 

addition to [H ], [Mg ] has also been included as a variable, especially because this ion may interfere with 
ATP species and so may influence J En through a change in [MgATP 2- ], which in turn would alter [CB]. 

Changes in [H + ] in the sarcosol are brought about mainly by two different mechanisms, which are 
both related to metabolic activity. One source of protons is manifest when metabolism is switched 

from rest to high power output. Fluxes in ATP consumption and production, J c ™ p and J APP , 

respectively, must then both increase to the same extent to reach a new steady state. During the 
adjustment, a phase of disturbed steady state occurs, during which both fluxes do not match. When 

power output increases, J C ° PP always leads J APP , i.e., there is an uncompensated ATP splitting until a 

new steady state is reached, at which point ATP production again equals ATP consumption. 

According to Alberty [20], this reaction is associated with proton production in dependence of [H + ] 
and [Mg 2+ ] (see (A6) and (A7)) for derivation of [H + ] changes and pH buffering). In addition, the CK 
and adenylate kinase (AK) reactions are involved, because these equilibria are also changed under 
these conditions and, as with ATP splitting, H + and Mg 2+ binding species are involved. Buffering of 
both ion concentrations is brought about mainly by sites intrinsic to the sarcosol. For Mg binding 
sites, an additional release of Mg 2+ by interfering [H + ] has to be expected. 
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Figure 5 A shows the time courses of rates of [H + ] changes. Interestingly, [H + ] production by ATP 
splitting is practically compensated by [H + ] consumption by the CK reaction. The contribution by the 
AK reaction is negligible. A similar behavior is found for Mg 2+ (Figure 5B). A concentration increase 
in this ion is mainly brought about by acidification. 

Figure 5. Time courses of [H + ] and [Mg 2+ ] during extreme power output. A: [H + ] fluxes; 
(brown) mainly LDH reaction and lactate transport; (red) ATP splitting; (blue) Jak, 
(yellow) Jck, (black line) resultant [H + ] flux; B: [Mg 2+ ] fluxes of the same reactions. 
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A second source of protons is given by the disturbance of lactate production by glycogenolysis (or 
glycolysis) and lactate efflux via lactate/H symport at the sarcolemma. Especially when lactate and H + 
accumulate in the glycocalyx (the outer aspect of the sarcolemma), the concentrations of these 
compounds also increase drastically in the sarcosol. This seems to be the main mechanism of 
sarcosolic acidification. 

Muscular fatigue at the cellular level can be defined as a phase of markedly reduced contractile 
performance, which largely recovers after a period of rest [38]. Because metabolites like creatine, 
ADP, Pi, H + , and lactate accumulate during conditions of fatigue in a similar way as can be observed 
during ischemia or hypoxia, which are known to be the result of impaired ATP production, it seems 
justified to suggest that the preconditioning for fatigue may also be initiated by a deterioration of the 
energy metabolism of the muscle fibers. Whenever ATP delivery does not match ATP consumption, 
such a situation may arise. 

These effects can be easily demonstrated with a simulation of glycogenolytic or glycolytic ATP 
production in the absence of mitochondrial metabolism (SJMoLYgen, see (A16)), which is related to the 
energy metabolism of fast muscle fibers. At 1.08 uM [Ca 2+ ] and a load of -1.5 x 10 4 J (constant 
glycogen content and glucose concentration [Glu] = 4.0 mM), efficiency of glycogenolytic ATP 

production is rj GLYgen = 0.722, that of glycolytic ATP rj GLY = 0.525. The higher efficiency is mainly 

caused by the stoichiometric coefficients of coupled ATP production of 3.0 and 2.0 for the 
glycogenolytic and glycolytic pathways, respectively. 

Under these conditions of high power output, metabolite concentrations change only moderately 
compared to resting conditions (at 1.06 uM [Ca 2+ ] and a load potential of -1.5 x 10 4 J/mol, [ADP] = 
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113, [Pi] = 8.32 x 10 3 , phosphocreatine concentration [PCr] = 9.7 x 10 3 , lactate concentration [Lac] = 3.0 x 
10 3 , [Mg 2+ ] = 832, and pH = 7.09). 

However, when a back pressure on glycogenolysis (or glycolysis) is produced by accumulated 
extracellular [Lac] e and [H + ] e , the flux through this pathway may become reduced. In addition, 
efficiency has been reduced by switching from glycogenolysis to glycolysis. The power output of ATP 
production is markedly reduced by these combined effects. As a result, the power of ATP production 
begins to fall, so that ATP consumption may overcome ATP production. Steady state cycling through 
ATP consuming and producing pathways can now no longer be maintained. 

Figure 6 shows that a first phase of slowly falling [MgATP 2- ] is followed by a phase of 
continuously enhanced reduction of this ATP species to low values ([MgATP ] = 230.0 uM; 
[PCr] = 1.6 uM). Immediately after reaching a minimum, a rapid recovery of [ATP] (up to starting 

2_|_ 

values) begins. [Mg ] shows a corresponding behavior. During the first phase it increases because of 
acidification, and then a sharp peak is produced by the onset of an extreme uncompensated ATP 
splitting (Figure 6). An increased [Mg 2+ ] may counteract the switch off of cross-bridge cycling and 
may aid recovery by increasing [MgATP "]. 

2 — 2~t~ 

Figure 6. Time courses of [MgATP ] and [Mg ] during development of fatigue, (red 
line) [MgATP 2- ]; (green points) [Mg 2+ ]. 
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Other parameters such as [PCr], [Pi], [Lac], and pH only partially recover under these conditions of 
extreme power output. An almost complete recovery, however, is possible under conditions of 
markedly reduced power output near resting [Ca 2+ ]. 

How this switch back to normal [ATP] is brought about can be seen from Figure 7. Not only have 

fluxes of ATP consumption and production, J C ™ P and J Pr ° p , become different now ( J € ™ p > J p ™ p ; 
Figure 7A), both fluxes of the cross-bridge cycle, J En and J Str , have also changed. These fluxes 

determine concentrations of [MH En ] and [CB], respectively. An increase in J En and a decrease in J Str 

would lower [MH En ] (whereby [CB] would be increased). Both concentrations always change 
reciprocally (Figure 7B). A P Str and A^ r are also affected. A P Str in particular is rapidly reduced until it is 



Metabolites 2012, 2 



687 



equal to - A Slr . Now all fluxes of the cycle must vanish, because the driving force of J Stl . has become 

zero. As a result, ATP consumption by cross-bridge cycling is switched off. 

Figure 7. Time courses of J c A ° T n p and J^™ , of [CB] and [MH En ], and potentials of the cross- 
bridge cycle during fatigue development. (A) (red) J A ° TP ; (blue) J P ™ P ; (B) (black) [CB]; 
(blue) [MH En ]; (C) (red points) A ATP ; (brown line) A P Str \ (blue line) A L s d tr ; notice that after 
7.5 s A P Str and A L s d tr become opposite and equal. 
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Under these conditions all myosin heads form cross-bridges, which however are unable to perform 
the power stroke, since the input force is equal to the opposed load force. In such a situation myosin 
heads may be bound to actin and may have dissociated H 2 P04~ and MgADP" similar to an isometric 
contraction, but in contrast to those latter conditions, equilibrium of forces is now brought about at a 

much lower load force (A P Str = ~ A£ = 0.375 x 10 4 J/mol at 1.08 uM [Ca 2+ ]). A load-dependent 

actomyosin splitting by MgATP 2- at the beginning of the stroke, that is uncoupling, is impossible 
under these conditions. So cross-bridge cycling with concomitant ATP consumption may be 
completely prevented. [ATP], therefore, can recover rapidly, even if the conditions leading to fatigue 
first remain unchanged. 

By this mechanism the fatigued skeletal muscle fiber is capable of protecting itself from the 
dangerous risk of irreversible cell damage. This seems to be necessary, since this cell type is 
voluntarily controlled without any protecting mechanism against an unbridled consumption of ATP, as 
is known to occur with other ATP coupled reactions such as, for instance, ion pumps. These are 
controlled mainly by the ion concentration that they are transporting. For example, when [Na + ] in the 
sarcosol is lowered by the Na/K pump to values below 10.0 mM, the reaction rate of this transport 
process is increasingly deactivated by the decreasing [Na + ], so that ATP consumption also is reduced. 

Such a protective mechanism is not known, however, for the cross-bridge cycle. Contractions with 
concomitant ATP splitting would be incessantly initiated, as long as firing of nervous impulses 
persisted. The voluntary muscle fiber would obey this parent command up to exhaustion or even up to 
cell death, if the fatigue producing mechanism were absent. Obviously it represents that special control 
mechanism which is necessary to protect voluntary muscle from dangerous ATP depletion during 
phases of high energetic demands. 
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The above results are obtained from a simulation, in which ATP production is confined solely to 
glycolysis. A whole muscle, however, is constructed from many types of functionally different fibers, 
with slow fibers having densely packed mitochondria, and fast fibers in which mitochondrial density 
can be very low. It is a known fact that especially fast fibers with a very low content of mitochondria 
and, therefore, a mainly anaerobic ATP production, are much more liable to be affected by 
overpowering than slow fibers. This may be brought about primarily by the preconditioning effects 
especially associated with this fiber type. At very high energetic demands, fast fibers produce much 
more lactate and protons through glycolysis than slow fibers, which can oxidise pyruvate by 
mitochondria. That is, the glycolytic ATP production rate of fast fibers may be decelerated by a back 
pressure, which may be generated by accumulating lactate and protons during high power output. In 
slow fibers with a high rate of oxidative glucose metabolism, such a back pressure cannot be produced 
as easily. Therefore, the metabolic changes leading to fatigue are simulated here with respect to fast 
fibers with a low resistance to fatigue. This weakness may be best demonstrated with an extreme fiber 
type, which can produce ATP solely by GLY. However, fibers completely devoid of mitochondria may 
not exist in vertebrate muscle. The results of this fatigue model, therefore, can only be taken as an 
approximation of real fast muscle fibers. Muscular contraction on the basis of NET has been treated 
theoretically by Caplan and Essig [13]. These authors explained the curvature of Hill's equation by an 
uncoupling. They defined the degree of coupling by q = \j yJ\ + 6 with O = a/F 0 = b/v max (6 between 

0.2 and 0.3). This latter equality is also fulfilled by the present model (0 = 0.309). However, in 
contrast to the approach of the above authors, the hyperbolic form of Hill's equation is produced here 
by introducing a Michaelis-Menten-like inhibition factor into the respective conductance, as already 
mentioned above. Moreover, uncoupling in the present model is produced through a load-dependent 
stroke shortening, which generates the maximum obtained with power plots. 

3. Methods 

Here the energy metabolism of a fast twitch muscle fiber is treated. That is, ATP production by this 
fiber type is solely brought about by metabolism of glycogen and/or glucose. Mitochondria are absent. 
Glycolytically produced [NAD re d] has to be reoxidised by the lactate dehydrogenase (LDH) reaction, 
and the lactate plus proton formed thereby is released to the extracellular space via Lac/H symport. 
Electrophysiological reactions at the cell membrane (sarcolemma) are omitted. Also, most reactions of 

2_|_ 

the sarcoplasmatic reticulum (SR) are not addressed. Only Ca pumping by the sarco/endoplasmatic 

2_|_ 

reticulum Ca ZT ATPase (SERCA) as an ATP consuming reaction is included in simulations besides 
several other reactions of ADP production (see SlMoLYgen (A16)) taken over from reference [1]. 

2+ 

Therefore, [Ca ] is treated as an adjustable constant. 

To determine the fractional fiber volume Vcai, a cylindrical geometry of the muscle cell is assumed. 
With radius R Ce ii = 25.76 urn, and a length L = 10 3 urn (fraction of whole fiber length), V Ce ii = 2.0847 x 10 6 
(am 3 or 2.0847 nL, and A Ce ii = 2.0847 x 10 3 urn 2 . From data of Aliev et al. [39] for the heart, the 
volume of the sarcosol, V c , can be determined by adding the mitochondrial to the fibrillar volume, 
yielding V c IV Ce ii = (321 + 195 + 55)/758.5 = 0.7528 or V c « 1.57 nL. Then a c can be obtained using 

a c = 10" 12 /(FxF c ) = 6.6024xl0" 9 uM/C (F = Faraday's constant, C = Coulomb). That is, to yield the 
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corresponding flux in uM/ms from an electric current entering the sarcosol, this current in 

in 

fA (= pSxmV =10 C/ms; pS = pico Siemens) has to be multiplied by a c . 

For calculation of force and velocity, the dimensions of the force generating cross-sectional area 
and the number of half-sarcomeres (HS) of the fibrils must be known. The contractile machinery of 
skeletal muscle fibers is organised in fibrils having diameters between 1.0 and 2.0 um, which are built 
up from in series sarcomeres connected by Z-discs over the whole length of a fibril, i.e., from end to 
end of the fiber. The functional unit is given by the HS. The principal filaments of an HS are the thick 
myosin and the thin actin filaments. In cross-sections, myosin filaments show hexagonal geometry. 
From this symmetry the fibrillar volume VFibr can be obtained. One hexagon is composed of six 
equilateral triangles of side length l Tr i = 41.0 nm [12] and equal angles of 60°. The area of a hexagonal 
fibril (or HS) of radius R Fibr = 18.0 x 41.0 = 738 nm is given by: 

(T8 / ) 2 

A Fibr = 6 - — ^-xsin (60°) = 1.415027 ^im 2 , and (19a) 
V Fihr = A Fihr xL Fihr = 1.415027 xlOVm 3 (=1.415027xl0~ 15 m 3 or 1415.027 pL) (19b) 
The total volume of fibrils is given by 0.866x Vceii (see reference [39]). The number NF&r is then given by: 

A^=^^ = 1.275842xl0 3 ( 19c ) 

For the determination of the total number of myosin heads of an HS, the number of myosin filaments 
of an HS must be known. The above hexagonal area of an HS can be constructed from equilateral 
triangles of lm = 41 .0 nm. A large triangle of the hexagonal HS with n-fold side length contains: 

MF m = 1 + 2 + 3- + n + (n + l) (20a) 

corner points representing myosin filaments. All points of the HS hexagon are given by: 

MF Hex = 6MF Tn -4 (» + 1) - (2 (n + 1) - 1) = 6MF Tn -6n-5 (20 b) 

For an 18-fold increase of lj ri from 41.0 nm to 738 nm (n = 18), the resulting HS hexagon contains, 
with MF T ri = 190, MF Hex = 1027 myosin filaments per fibril. 

Myosin filaments contain 294 myosin heads per half filament. In an HS, thus 294 x MF Hex , and in 
the whole fiber: 

N^ = 294xMF hex xN Fibr = 3.852252xl0 8 (20c) 

myosin heads are contained. 

Their concentration is obtained by dividing by Na (= 6.022142 x 10 particles per mol), and by the 
volume of the water diffusible space in the filament lattice containing myosin heads of all HSs of the 
cross-sectional area of a fiber, Vlm- This volume is given by: 

Kl" = f WDS (A Flhr x f MH l HS x N Flhr ) = 972.453nm 3 or 0.972453 pL (20 d) 

(Ihs = HS length at rest =1.1 \im,fMH = length fraction of myosin filament containing myosin heads in 
terms of HS ([12]) = 0.62364, f WDS = volume fraction of water-diffusible space in the filament lattice 
volume = 0.7852). 

The myosin head concentration is given by: 
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N tot 

[ MH L = TT-^«= 657.8 uM. ( 20e) 

At steady state, the concentration of stroking cross-bridges, [CB] = [MH] M - [MH] ([MH] = the 
remaining myosin head concentration given by the simulation) is determined by two fluxes (see 
Results). It is adjusted to about 25% of [MH] tot at [Ca 2+ ] = 1.08 uM and a load potential A L = 3.0x1 0 4 



J/mol. The force generated by [CB] is given by: 

]CB] 
[MH] 



F = F c B ^ Hjt - t - ( 21a) 



F C b is the force of one single cross-bridge. It is obtained here from the stroke input potential A P (see 
Results) under these conditions, by assuming that the stroke length ls tr = 12.0 nm. Then Fcb is given by: 

A p 5.819697 xlO 4 onc ^, MA ^ 

per cross-bridge, yielding F = 0.7735 mN per fiber, and a tension of 3. 71 x 10 5 Pa (related to A Ce ii, 
Pa = Pasqual). The force is generated at 1.08 uM [Ca 2+ ] by about 25% of myosin heads ([CB] = 164 
uM). All force-generating cross-bridges are contained in the parallel HSs of the cross-sectional area of 
1 027 fibrils of one muscle fiber. 

Simulations were solved using Mathcad® 14.0 or 15.0 M01 1 solver AdamsBDF. The programs were 
run under Microsoft® Windows 7 and XP Professional. 



4. Conclusions 



Cycling between coupled reactions occurs, especially in energy metabolism. It is shown that the 
overall resistance of such cycles must vanish, and that the resistance or conductance associated with 
the negative output affinity of a coupled reaction also has to be negative. The following may be 
illuminating: The entropy change of a spontaneously proceeding reaction (4A ) i s always positive. 
When a reaction is forced against spontaneity, A r S i must become negative. All reaction parameters 

associated with entropy changes, like affinities and conductances, must inevitably follow a sign change 
of A r S i whenever such a change occurs. This is not a contradiction to Ohm's law, but a consequence 

of the phenomenological definition of a conductance through +L = J/+A . It can be concluded that the 
occurrence of negative conductances is realised not only with biochemical reactions in living cells, but 
represents a fundamental concept of coupled processes. 

This concept is realised here for the cross-bridge cycle. The reactions of the cycle are described on a 
thermodynamic basis using the kinetic approach of enzyme-catalysed reactions. Hill's equation for 
muscular performance can be derived on this basis. However, uncoupling has to be introduced to yield 
a maximal efficiency of power output. Here the uncoupling mechanism is not an accidental process 
during energy transduction, but a necessary interference during force generation, which ultimately 
produces an isometric contraction. 

Although mechanical acceleration may also be possible on a cellular basis by changes in sarcosolic 
[Ca 2+ ], it seems highly unlikely, however, that this may be sufficient to allow normal locomotion of a 
subject. Only the control by the nervous system can bring about coordinated actions of several muscle 
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fibers, groups of fibers, or even several different muscles. In this way, accelerated and decelerated 
motion becomes possible. To achieve this, the number of force generating cross-bridges is varied by a 
change in cross-sectional area, that is, by altering the number of fibers recruited for contraction. 
Thereby the locomotion at high efficiency or maximal power output can be controlled by will. Also, 
isometric contractions are indispensible for coordinated actions. They are produced by reducing the 
cross-sectional area to such an extent that a load dependent uncoupling is initiated to stop 
fiber shortening. 

In many species nervous control of muscles is not a capability which is present from birth on. To 
reach a certain level of adroitness an individual has to learn — often during a long lasting phase of 
exercise — to control muscle action by coordination. 

Appendix 

Negative Resistances in Simple Electric Circuits 

Reactions occurring in a common car battery can be considered as coupled. A redox reaction is 
started by introducing a catalyst (the electrodes), which couples the affinity A R of the redox reaction to 
the formation of an electrical potential difference A<f> at the electrodes. Under open circuit conditions 
the reaction proceeds rapidly to equilibrium, at which A R +A(f> = 0. Taking A R as the positive input 
force, then A<j> must be negative. A R can be expressed in electrical units using E = A R /zF (E = 
electromotive force in Volt V, z = charge number, F = Faraday constant in Coulombs/Volt, A0 in V). 

The coupled flow of charges (electrical current in ampere A) is then given by: 

I = L c (A</> + E) ( A1) 

E and A(j) correspond to the input force A 2 and output force Aj, respectively. L c is the coupling 
conductance, and R t = \IL C represents the inner resistance of the battery (R's are given in Q., L's in 
1/H). The partial conductances L c i and L C 2 (see equations 2e and 2f) are given by: 

A0 + E A0 + E 

cl = 4 ~A0~ ' c2 = 4 ~~E~ (A2) 

In a simple electric circuit consisting of a battery (E = 12 V), an inner resistance R { (0.2 f2) and an 
outer resistance R e (0.4 f2), the current / is given by E - 1 + R e ) = 0 , yielding: 

/ = - = 20.0 A , and from (Al), A(j> = I ~ L - E = —I = -8 V . 
R+R L L 

i e c e 

The outgoing voltage delivered to the outer circuit is then given by U e = —A(j> = +8 V. 

L c i and L c2 are given by -5/2 and 5/3 f2 _1 , respectively, or Rn = -0.4, and Ra = 0.6 Q, which fulfils 
Ru + R i2 = Ri = 0.2 D.. Setting R e = 4.0 D. yields R a = -4.0, R i2 = 4.2 Q., and again R t = 0.2 Q.. 
Moreover, in the electric circuit the overall resistance Rn + R e vanishes. It should be noticed that 
partial resistances are not constant, although R t is a constant. They depend both on R e and L e , 
respectively. 

The total dissipation function of reactions in the battery and of the outer circuit is given by: 
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®drc = 0 1 +0 2 +& e = l{A<fi + E + U), 0 circ =IxE (240 J/s or 0.24 kW), or (A3a) 
®drc = hi x E 1 (l/0.6x 144 = 0.24 kW) (A3b) 

This latter result means that the total entropy production of the circuit is given by that of the redox 
reaction. This is also valid with respect to heat production. 
Power output as a function of R e is given by: 

p - F 1 

out (R e+Ri y < A4a ) 

Maximal power output is reached at R e = R t = 0.2 Q, P™f = 0.18 kW. Because P out can also be 

expressed as P out = -L cl xA(f> 2 , this yields withZ, = 1/0.2 = -(-5.0)( -6) 2 = 0.18 kW. 

When a second battery is added to the circuit in such a way that E 11 is directed against 
E 1 (E n = -10.9 V, Rf = 0.5 Q.), A<f> n is now positive and E n negative. The outer resistance R e stands 

for the resistance of the wires connecting both batteries. In this constellation, - Atp 1 is no longer equal 
to U e . Now - (At/) 1 + U e )= - A(J) Ie = A(J)" is valid. Consequently, R e also has to be added to R* 
yielding R] e = Rj + R e . 

/ = ^T^ =1M - U. = R.I=0AV (A4b) 

from: 

I = ±(Af + E 1 ) , and / = ±-(E a + Af ) (A4c) 

A(j> le = -1 1.4 V, and Af = 1 1.4 V is obtained. 
Partial conductances are given by: 

L' cl = L 1 : ^jf- = - 0.0877 and L" c2 = L» E " = 0.0877 I/O. (A4d) 

Again, the overall resistance of the electric path in the circuit is zero. The total resistance 
Bj +R e +R? = 1.1 £2, therefore, is also given by R^ + Rjj = 12 + (-10.9) = 1.1 £1. 

ATP, ADP, and Pi Species as Functions of [ft] and [Mg 2+ ] 

ATP species, including MgATP , are calculated according to the methods of Alberty [20]. When 
respective constants are known, which are dependent on temperature and ionic strength, so-called 
polynomials can be formulated, from which several parameters like species concentration, K' 
(biochemical equilibrium constant), or [H + ] and [Mg 2+ ] binding can be taken. 
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ATP splitting by myosin ATPase is formulated here for the species MgATP 2- , MgADP", and 
H 2 P0 4 ~, to make the reactions of the cross-bridge cycle directly dependent on these compounds. The 
reaction in chemical notation form is given by: 

MgATP 2 " + H 2 0 ^ MgADP" + H 2 P0 4 

The equilibrium constant for the above reaction is given by: 

K reQ = K refl x [ p ^3-l[^-] = 6 . 8 28 x 10 4 (A5a) 
P -P x ° at 

r MgATP2- A1PA - ~Mg + ~\ (A5b) 

P -P x ^ 

-*- \A rr Ar\D -*- Anon ^ i 



MgADP- 1 ADP3- * J" Mg2+ j (A5 C) 

K refl (= 6.267 x 10 5 ), K3at, K3ad, P A tp 4 -, and Padps- were taken from [1]. At given [H + ] and [Mg 2+ ] 
values, the corresponding K' (= 4.9687 x 10 5 , pH = 7.1, [Mg 2+ ] = 800 uM) is identical to formulations 
with other reference constants. 

[Ft] and [Mg 2+ J Buffering 

[H + ] buffering of SMFs is treated here analogously to VMs (see [1]). It is given by: 

4H+]_^_ (A6a) 
dt \ + /c BuZjii 

At very high ATP consumption, which may not be matched by ATP production differences in 
proton binding between substrates and products of reactions like ATP splitting, CK and AK reactions 
may also contribute to pH changes. According to [20], such changes for [H + ] and [Mg 2+ ], respectively, 
are given by: 

dlnK' , .- dlnK' 
A#h = r TT+ n > and AN Mg = 



dln[H + Y r Mg dln[Mg 2+ ] ( A6b ) 
- dlnP r , , dlnP 

N " = JifiF] ■ and N >* {HMg) = 5h[Sp] < A <*> 

as a result, 

= 7- BUT BU X (Z V i J i ~ AN H (J ATP + J C K + JaK )) ■ (A6d) 



-f— ^ xBU'h i s identical to k b h u if binding sites contain only a single site with only one proton 
d I II ' 



dissociation constant. 
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For [Mg 2+ ] buffering, it is suggested that during short time intervals Mg 2+ transport reactions across 
membranes can be neglected. Only intrinsic binding sites including [ATP] are present and, as with 
[H + ] changes, [Mg 2+ ] changes induced by ATP splitting, the CK reaction, and the AK reaction have 
been addressed. [Mg ] buffering can be expressed as: 

- L ^" J + 1 ~ = ~AN Mg (JZ + Jck (A7a) 

4^] hound _ dN Mg 4Mg 2+ ] 

dt i[Mg 2+ ] dt (A7b) 

jrMg 2+ l AN M , v 

_L & -I - r M g x { T con + T + T \ 

M rIAT V CK AK J /A „s 

dt g^_ x r BU p (A7c) 

j[Mg 2+ ] L jMg 

In addition, Mg 2+ binding depends on [H + ]. A decrease of pH can liberate magnesium ions from 
intrinsic binding sites and from the predominant ATP species MgATP 2- . The H + and Mg 2+ dissociation 
constants of both binding sites are set to the values of a simplified Patp4— The total concentration of 
Mg 2+ binding sites, [BU]™ , is adjusted to 9.0 mM plus a variable [ATP]. The change of [Mg 2+ ] is 

given then by: 

4Mg 2+ ]_ AN Mg 



( jcon j r \ 

y ATP CK AK J 



J L M § J (A7d) 



-^r([BUE + [ATP])xlM. 



In simulations, instead of complete d[H + ]/dt, only those fluxes producing or consuming protons are 
considered, because changes of [H + ] depend mainly on these fluxes (see Figure 5 A). 

[Mg ] is introduced as a variable only in those simulations that deal with muscular fatigue. 
Because changes of [Mg 2+ ] depend mainly on acidification, and pH does not change markedly even 

under conditions of high power output, this variable is set constant to 800 uM for all other simulations. 

In the above equations, methods of calculus are used so formulas can be held compact. In 
simulations, however, these equations must be incorporated in an explicit form, which often results in 
very voluminous expressions. 

Simulation of Glycogenolysis and Glycolysis 

Most flux equations of glycogenolysis are congruent with those of a simulation of glycolysis given 
in [1]; they are taken over from that article. Glucose-6-phosphate (G6P) formation by hexokinase (HK) 
and glycogen phosphorylase is now included. The new flux equations used here are as follows. 
Flux through glycogen phosphorylase: 
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Phosph 



[AMP] 2 

(C S "*) 2 +[AMP] 



2 Phosph 



x RT In 



[Pi] 
[G6P] 



k: 



Phosph 



^w/, =4xlO- 3 (^M/ms)x(mol/J), K?* = 2.0 liM, K' = 0.286; 



glucose - 6 -phosphate isomerase, 

[G6P] 2 



J GPL 



(KT) + [G6P] 



[F6P] ' 



GPl 



-'GPl 



2xl0" 2 (uM/ms)x(mol/J), = 300 uM, K^ P1 = 0.276; 



lactate dehydrogenase, 



J ' LDH 



[Pyrf 



K*)+[Pyr] 



xRT/n 



[Lac] [NAD ox ] 



LDH 



j max 
LDH 



2.4xlO" z (uM/ms)x(mol/J), =50 uM, K' LDH =2.497x10 



lactate/proton cotransport, 



J 



[Lac] 



RT 

+ [Lac] iac F 



[Lac] c [H + ] c 
[Lac] e [ H +] 

G™ x = 2.866x 10 8 pS (pico Siemens = 10" 8 ^ -1 ), K^ ac = 17 mM; 
Na + /H + exchange, 

KIM/ 

K1[H*]J 



^Ata// _ a c 



\ + e 



1 RT 

h 0 3-[h + P G ™ hX l " 



s H 



= 10 5 pS, H 05 = 0.1 uM, S [H + ] = 0.004 uM; 



anion exchange reaction, 
1 



J AnEx a c 



[HC0 3 ] 



1 + e 



*M (K-f +[HC0 3 ] 



2 / 
RT 

1 max 7 

2 <^ n£x x — In 



[Cl] e [HC0 3 ] c 



[Cl] c [HCO3] 



e / 



G;„" x = 10 4 pS, H 05 = 0.05 uM, S 05 = 0.008 uM, K™ ex = 13.0 mM 



The energising flux of the cross-bridge cycle is given by: 

f ( 



J En 



[CB] 
[CB 0 ] 



1+e 



Zmax 
c„ X 



A L05 J coir A Slr 



A„+RTln 



[ATP] ] 

[CB] X (p ) 



IS-ref 

K3at 



[Mg 2+ ] 
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L m E a n x = 6.138xl0" 2 (uM/ms)x(mol/J),.£ w = ([CB t ]-[CB 0 ])/([CB t ]-[CB]), [CB t ] = 656 uM, [CB 0 ] = 



492 uM, e = 24.0, A L0 5 = 3.0xl0 4 (J/mol), S L = 2.0x1 0 4 (J/mol), , A En 



-5.866729xl0 4 



Kg f = 5.2 x 10 8 



0.326/ 



f[Ca 2+ ]-0.55^ ^ 



l + e 



+ 1 

J J 



The factor [CB]/[CB 0 ] corrects L m E a * for changes of [CB]; the second factor is introduced to damp 



changes of [CB]. K r / is not constant, but depends on [Ca 2+ ]. 
The stroke generating fluxes are given by: 



_ [CB t ]-[CB] 

J Str —' 



1 



K 



CB 



[CB t ]-[CB 0 ] 



l + e 



Ca,,,- Ca 



^1 + Jcorr A Str 



Ld Str 



f f 



A 



Ld 
'Str I 



l + e 



+ ■ 



A 



Ld 
'Strl 



+ 1 



Ld 



X fcorr A Str 



+ RT/n 



l + e v 12 
[CB][c°] iW^ 4 -K3ad 



[ADP][Pi] [Mg 2+ ] 



K ref 



_ [CB t ]-[CB] 



[CB t ]-[CB 0 ] 



l + e 



Ca 2 



^7 + fcorr ^Str 



Ld Str 



f 

f A Ld + 

Jcorr ^Str T 



/I,, 



/cog; ^flr rt 0 
Sii 



+ ■ 



A 



Sir 2 



X 

V v 



RT In 



l + e v " ; l + e v 
[CB][c°] P ADP3 _P HiPOA _K3 a d 



+ 1 



[ADP][Pi] [Mg 2+ ] 



K 



ref 



+ 



(A15a) 



(A 15b) 



Both fluxes are identical as long as uncoupling is absent. The first factor corrects Z™ x = 4.6 x 10 



(uM/ms)x(mol/J) for changes of [CB t ] - [CB]. The second factor introduces [Ca 2+ ] dependence of L Str 



2_j_ 

The third factor is responsible for the hyperbolic character of the flux equation at constant [Ca ] with 



Kf B = -1.8 x 10 J/mol, which represents the inhibition constant, K r / = 1.310889 x 10 4 . X values are 

not independent; this interdependency is given in Results. Uncoupling is formulated to occur in two 
steps, expressed by A F Strl = 0.15 and A p Strl = 0.85. 

SlM GLYgen (A16) 



d [ ADP] 



dt 



J\Vser + ^Wb + ^ Str + ^ GK + ^ PFK ^ Ga ^ AK ^ CK • 
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_L J = / + /+/_/_/ 

^ ^ ^ser Wb Str ° Ga J Phosph 



!&1 = J + J -J 

^ HK Phosph J GPI ■ 



J[F6P] 



dt 



^GPl J pfk ■ 



4FBP] 



dt 



J ' PFK ^ ' Al 



d J2^1 = j + j -j 

dt M T > Ga 



^[DHAP] 



dt 



~ J Al ^Ti 



4Pyr] 



dt 



J Ga JlDH 



d[Lac] 



dt 



J ' LDH J Lac 



4NAD red ]_ 



dt 



J T Ga J J Lac • 



4coj_ 



dt 



JcA JcOlc ■ 



</[hco 3 ] 



dt 



' JcA J AnEx 



d[CB] 

dt 

j[h + ] 



_ J En Jstr ■ 



dt 1 + 



J[ADP] 



dt 
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4Mg 2+ ] = _AN Mg J djAm] 



dt 1+ f Nilg -. x [BU] tot 

^[Mg 2+ ] L J 



V 



dt 



Mg 



+ 



-^([BUE + [ATP])xlM. 



The above set of differential equations without a variable [Mg 2+ ] ([Mg 2+ ] = 0.8 mM = const.) is 
used to calculate the various points of figures (Figures IB, 2, 3, and 4) for a given [Ca ] and various 
loads. As already mentioned, [Mg ] is introduced as a variable only for conditions of very high power 
output leading to fatigue. From the output of the simulation many more variables, as shown here, can 
be obtained as functions of time, which may often be helpful in understanding underlying mechanisms. 
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